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^ The ultrahigh energy neutrino cross section is a crucial ingredient in the calculation of 

the event rate in high energy neutrino telescopes. Currently there are several approaches 
which predict different behaviors for its magnitude for ultrahigh energies. In this con- 
tribution is presented a summary of current predictions based on the non-linear QCD 
^ ' evolution equations, the so-called perturbative saturation physics. In particular, pre- 

, dictions are shown based on the parton saturation approaches and the consequences of 

■ geometric scaling property at high energies are discussed. The scaling property allows an 

' analytical computation of the neutrino scattering on nucleon/ nucleus at high energies, 

' providing a theoretical parameterization. 

1. Introduction 

The investigation of ultrahigh energy (UHE) cosmic neutrinos provides an oppor- 
tunity for study particle physics beyond the reach of the LHC As an example, 
nowadays the Pierre Auger Observatory is sensitive to neutrinos of energy > 10^ 
\^ • GeV^. A crucial ingredient in the calculation of attenuation of neutrinos traversing 

the Earth and the event rate in high energy neutrino telescopes is the high energy 
neutrino-nucleon cross section, which provides a probe of Quantum Chromodynam- 
ics (QCD) in the kinematic region of very small values of Bjorken-x. The typical x 
value probed is a; w m^^ /2mNE^, which implies that for w 10^ — 10^° GeV one 
have X w 10"^ - 10"^ at ~ GeV^. This kinematical range was not explored 
by the HERA measurements of the structure functions 

The description of QCD dynamics in such very high energy limit is a subject of 
intense debate SI. Theoretically, at high energies (small Bjorken-a;) one expects the 
transition of the regime described by the linear dynamics, where only the parton 
emissions are considered, to a new regime where the physical process of recombi- 
nation of partons becomes important in the parton cascade and the evolution is 
given by a non-linear evolution equation. This regime is characterized by the lim- 
itation on the maximum phase-space parton density that can be reached in the 
hadron wavefunction (parton saturation), with the transition being specified by a 
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typical scale, which is energy dependent and is called saturation scale Qs- More- 
over, the growth of the parton distribution is expected to saturate, forming a Color 
Glass Condensate (CGC), whose evolution with energy is described by an infinite 
hierarchy of coupled equations for the correlators of Wilson lines SI. In the mean 
field approximation, the first equation of this hierarchy decouples and boils down 
to a single non-linear integro-differential equation: the Balitsky-Kovchegov (BK) 
equation . Experimentally, possible signals of parton saturation have already 
been observed both in ep deep inelastic scattering at HERA and in deuteron-gold 
collisions at RHIC 

Currently, there are predictions of the neutrino nucleon cross sections with struc- 
ture functions constrained by HERA data are based on linear dynamics using 
DGLAP or an unified DGLAP/BFKL evolution, or phenomenological models that 
resembles the expected behavior predicted by the non- linear QCD dynamics!^ (i.e., 
the proton structure function saturating the Froissart bound at asymptotic energies, 
F2 cx ln^(l/a;)). As a general feature, the nonlinear QCD dynamics predicts sizable 
suppression of UHE neutrino cross section in comparison with st andar d approaches. 
Here, we summarize the main results of works presented in Refs. ECni j^gf. ^the 
contribution of non-linear effects was estimated considering distinct phenomeno- 
logical models based on satu rati on physics. An update on those calculations have 
been done recently . In Ref . the geometric scaling property (which is a natural 
consequence of the asymptotic solutions of the nonlinear QCD evolution equations) 
is considered to obtain an analytical parameterization for the UHE neutrino cross 
sections. In what follows we introduce the theoretical and phenomenological tools 
and present the main results and predictions 



2. The UHE neutrino cross section and nonlinear perturbative 
QCD approaches 

Deep inelastic neutrino scattering is described in terms of charged current (CC) and 
neutral current (NC) interactions, which proceed through and Z° exchanges, 
respectively. The total cross sections are given by: 

2 



CC. NC 
''uN 



jQ2/gXS TT 
2 



Mf 



1 + (1 - y) pCC, NC _ V_pCC, NC 



(1) 



where En is the neutrino energy, s = 2MEu with M the nucleon mass, y — Q"^ /(xs) 
and Q^jjjj is the minimum value of which is introduced in order to stay in the 
deep inelastic region. In what follows we assume Q^in = 1 GeV'^. Moreover, G_f is 
the Fermi constant and Mi denotes the mass of the charged of neutral gauge boson. 

In the QCD improved parton model the structure functions Fi{x,Q^) are ex- 
pressed in terms of the parton distributions on the nucleon, which satisfy the 
DGLAP linear dynamics. On the other hand, an efficient way of introducing 
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Geometric scaling in ultrahigh energy neutrinos and nonlinear perturbative QCD 3 

non-linear effects is the description of the structure functions considering the color 
dipole approach in which the DIS to low x can be viewed as a result of the interac- 
tion of a color qq dipole which the gauge boson fluctuates In this approach the 
F2 ' structure function is expressed in terms of the transverse and longitudinal 
structure functions, = + which are given by 

F^^f'^^{x,Q')^^j\z J dV|<f(r,z,Q2)|V^^(r,x) (2) 

where r denotes the transverse size of the dipole, z is the longitudinal momentum 
fraction carried by a quark and ^^'^ are proportional to the wave functions of 
the virtual charged of neutral gauge bosons corresponding to their transverse or 
longitudinal polarizations. Explicit expressions for ^'^'^ are given, e.g., in Ref. 
Furthermore, tr''*^ describes the interaction of the color dipole with the target. In 
next section we will discuss some models for a'^^^, based on the non-linear QCD 
dynamics, which describe the current HERA data. As a comment, the DESY-HERA 
measurements of the structure function at low - x (a; w 10~^) are for very low values 
of (Q^ ^ 1 GeV^), which implies that small x extrapolations of the parton 
distributions are necessary to estimate a^^N above E^^ w 10^ GeV. 

As mentioned above, UHE neutrino-nucleon cross section accesses very large 
values of and very small values of Bjorken x. The usual linear perturbative 
QCD evolution equations predictions are based on the linear DGLAP and/or BFKL 
equations which implies a power increase with the energy of the neutrino-nucleon 
cross section that eventually violate the Froissart bound. A easy way to introduce 
consistently the unitarity constraints on the theoretical estimate of the ultrahigh 
energy behavior of the neutrino-nucleon cross section is to express the structure 
functions in the dipole approach [Eq. ([2])] and to consider the state-of-art of the non- 
linear QCD dynamics: the Color Glass Condensate formalism. In this formalism, 
the dipole - target cross section a'^'^P can be computed in the eikonal approximation, 
resulting 

a'^''P{x,r) = 2 J dHj\fix,r,b) , (3) 

where Af is the dipole-target forward scattering amplitude for a given impact pa- 
rameter b which encodes all the information about the hadronic scattering, and 
thus about the non-linear and quantum effects in the hadron wave function. It is 
useful to assume that the impact parameter dependence of TV can be factorized 
as M{x,r,b) = J^{x,r)S{b), so that cr''*P(a;, r) = aoAf{x,r), with ctq being a free 
parameter related to the non-perturbative QCD physics. The Balitsky-JIMWLK 
hierarchy describes the energy evolution of the dipole-target scattering amplitude 
Af{x, r). In the mean field approximation, the first equation o f thi s hierarchy decou- 
ples and boils down to the Balitsky-Kovchegov (BK) equation^EI xhe dipole-target 
cross section can also be calculated considering phenomenological parameterizations 
for J\f{x, r) based on saturation physics, which provide an economical description of 
a wide range of data with a few parameters. This approach has been considered in 
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Refs. ISCH In general, the dipole scattering amplitude is modeled in the coordinate 
space in terms of a simple Glauber-like formula as follows 

1 



Af{x, r) = 1 — exp 



(4) 



where 7 is the anomalous dimension of the target gluon distribution. The main 
difference among the distinct phenomenological models comes from the predicted 
behavior for the anomalous dimension, which determines the transition from the 
non- linear to the extended geometric scaling regimes, as well as from the extended 
geometric scaling to the DGLAP regime SI. 

One important property of the nonlinear perturbative QCD approaches for high 
energy deep inelastic ep{A) scattering is the prediction of the geometric scaling. 
Namely, the total j*p{A) cross section at large energies is not a function of the 
two independent variables x and , but is rather a function of the single variable 
TA — Q^/Qs^^ A- usual, is the photon virtuality and x the Bjorken variable. 
Geometric scaling is the exact asymptotic solution of a general class of nonlinear 
evolution equations and it appears as a universal property of these kind of equations. 
The saturation momentum QL,a(2;; ^) ^ A"' x'^ {a ~ 1/3, A ~ 

0.3) is connected with the phenomenon of gluon saturation. In principle, geometric 
scaling is predicted to be present only on process dominated by low momenta. 
However, it is known that the geometric scaling is preserved by the QCD evolution 
up to relatively large virtualities within the kinematical window Qsat(^) < < 
Qsnti^) / -^QCD- That is, the scaling property extends towards very large virtualities 
provided one stays in low-x. This kinematical window is further enlarged due to the 
nuclear enhancement of the saturation scale. These facts have direct consequences 
in the behavior of UHE neutrino cross section. 

As similar scaling holds for the massive boson-nucleon cross section, we are able 
to compute analytically the neutrino cross sections. This feature was verified to be 
true in Ref. ^1^, where geometric scaling was shown to occur in the small-x charged 
current neutrino data. In this sense, we are able to construct a theoretical prediction 
which is model independent. They are as follows 

' , (5) 



where Af[i) are overall normalizations, 2 Eire numerical constants with 
CC, NC, Wscai — bX and a — b/S. This implies in a mild power-like rise Wscai — 0.2 
for the neutrino cross section in contrast with other theoretical approaches. The 
nuclear dependence is approximately linear, (J^^^^°^ oc A(t"'^'^'®°", once 6 ~ 5 and 
hence a ~ I. The remaining constants are given by (see Ref. for details). 



aoGl M2 \ /ax"— \ B 



StT^X j \ b J ttcm J2f^f 



Cf^ =7r(2m^)'^— (M^_2)"'=^^'csc(^i.,eai) (1 - i^scai) , C^'^ = tt {M^^z)'" CSC (nb) [l - b) . 
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Fig. 1. Comparison among perturbative QCD approaches: unified DGLAP-BFKL, NLO DGLAP 
and a phenomonological model based on the Color Glass Condensate physics (CGC). The analyt- 
ical geometric scaling parameterization is also presented. 

where one uses the notation Bqc — 4, -Bnc = ^chirai and j/gcai = b — Wscai- Numer- 
ically, this gives a total cross section cr*°^ = 1.48 x 10^^^ A°'{E^/GeVf '^^'' cm^. 
The cross section above can have imphcations for neutrino observatories because 
experiments are planned to detect UHE by observation of the nearly horizontal air 
showers in Earth coming from neutrino- air interactions ^1^. A reduced cross section 
produces a smaller event rate for such neutrino-induced showers and could com- 
promise the detection signal. However, the rate of up-going air showers initiated 
by muon and tau leptons produced in neutrino-nucleon reactions just below the 
surface would increase, being possibly larger than the horizontal air shower rate. 

3. Results and summary 

Let us present the comparison between the predictions of the linear approaches 
(NLO DGLAP and unified DGLAP/BFKL) and non-linear QCD approaches (a 
phenomenological model for the Color Glass Condensate and the geometric scal- 
ing prediction). In Fig. [1] the energy dependence of the neutrino nucleon GC cross 
section predicted by the linear and nonlinear perturbative QCD approaches are 
compared. As expected from the solution of the DGLAP equation at small- a;, the 
NLO DGLAP (dot-dashed curve) and unified DGLAP/BFKL (long-dashed curve) 
results predict a strong increase of the cross section at ultrahigh energies. Such es- 
timates are largely used to estimate the event rates in neutrinos telescopes. Clearly, 
the CGC calculation predicts an increasingly suppressed cross section compared to 
its linear QCD counterpart. 

In Fig. [T] the solid line represents the result of geometric scaling parameteri- 
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zation, Eq. ([S|). The suppression in the cross section is even higher than for the 
CGC phenomenological models reaching a factor two compared the the Unear QCD 
approaches. The present resuhs demonstrate that the determination of a^N can be 
useful to constrain the underlying QCD dynamics. In principle, this cross section 
could be constrained at high energies by studying the ratio betwee n qu asi-horizontal 
deeply penetrating air showers and Earth-skimming tau showers ^1^. 

As a summary, the detection of UHE neutrinos may shed light on the obser- 
vation of air showers events with energies in excess of 10^^ GeV, reveal aspects of 
new physics as well as of the QCD dynamics at high energies. One of the main 
ingredients for estimating event rates in neutrino telescopes and cosmic ray obser- 
vatories (e.g. AUGER) is the neutrino - nucleon cross section. We have examined to 
what extent the cross section is sensitive to the presence of new dynamical effects 
in the QCD evolution. We compare the predictions of several approaches based 
on different assumptions for the QCD dynamics. In particular, we have compared 
the phenomenological models based on nonlinear perturbative QCD with the usual 
NLO DGLAP evolution. Moreover, we compute a model independent cross section, 
using the experimental verification of the geometric scaling property. The results 
demonstrate that the current theoretical uncertainty for the neutrino-nucleon cross 
section reaches a factor two or higher for neutrino energies above 10^^ GeV. 
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